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ABSTRACT 



Context. Rotation Measure synthesis of the Westerbork Synthesis Radio Telescope (WSRT) observations at /I ~ 2 m of the FAN 
region at /=137°, b-+l° shows the morphology of structures in the ionized interstellar medium. 

Aims. We interpret the diffuse polarized synchrotron emission in terms of coherent structures in the interstellar medium and the 
properties of the interstellar magnetic field. 

Methods. We performed statistical analysis of the polarization data cube obtained through Rotation Measure synthesis. For the first 
time, cross-correlation is applied to identify and characterize polarized structures in Faraday depth space. Complementary information 
about the medium are derived from Hor emission, properties of nearby pulsars, and optical polarized starlight measurements. 
Results. We find an overall asymmetric Faraday dispersion function in a Faraday depth range of [-13, +5] rad m~^, which is peaked 
around -1 rad m~^. Three morphological patterns are recognized, showing structures on scales from degrees down to the beam size. 
The first structure is a nearby synchrotron emission component with low Faraday depth, filling the entire field of view. The second 
pattern is a circular polarization structure with enhanced (negative) Faraday depth, which has the same morphology as a low-emission 
region within the third component. This third component is interpreted as the background in which the circular structure is embedded. 
At low Faraday depth values, a low gradient across the imaged field is detected, almost aligned with the Galactic plane. Power spectra 
of polarized structures in Faraday depth space provide evidence of turbulence. 

Conclusions. A sign reversal in Faraday depth from the nearby component to the circular component indicates a reversal of the 
magnetic field component along the line of sight, from towards the observer and nearby to away from the observer at large distances. 
The distance to the nearby, extended component is estimated as < 100 pc, which suggests that this structure corresponds to the Local 
Bubble wall. For the circular component, various physical interpretations are discussed. The most likely explanation is that the circular 
component seems to be the presence of a nearby (~ 200 pc away) relic Stromgren sphere, associated with an old unidentified white 
dwarf star and expanding in a low-density environment. 
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1. Introduction 
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The study of the magneto-ionic properties of the interstellar 
medium (ISM) by multi-frequency pol arimetry started in 1960 
with the Dwingeloo radio telescope ( Westerhout et al. 1962|) . 
A breakthrough came much later when |Wieringa et al. ( [1993| ) 



. ^ acknowledged the importance of small angular scale structures 
K> (down to a few arcmins) in polarized intensity (PI) or polariza- 
^ tion angle with no detection in total intensity. Such a fine struc- 
ture was explained in terms of Faraday rotation of the diflTuse, 
linearly polarized, Galactic synchrotron radiation background by 
a highly structured foreground magneto-ionic medium. The cos- 
mic rays propagating through the Galactic magnetic field are re- 
sponsible for the radio synchrotron emission, and the magneto- 
ionic medium is primarily the warm ionized ISM. As a result 
polarimetric observations of the Galaxy provide insight into the 
ionized ISM and the Galactic magnetic field. 

The FAN region is a highly polarized, extended region in the 
second Galactic quadrant at low Galactic latitudes. Its polariza- 



t ion pr operties were first investigated by [ Bingham & Shakeshaft 
( jl967| ), who discovered a circular structure in the rotation mea- 
sure (RM) map. From higher angular resolution observations, 
Verschuur| ( |1969[ ) suggested a connection of the observed struc- 
ture with the B2IV:e star HIP 15520 that is located close to the 
centre of the circular structure. 



[ Haverkorn et al.| ( |2003| ) presents multi-frequency polarime- 
try radio observations at five wavelengths in the range of A = 
0.80 - 0.88 m. They point out a ring structure in PI along with an 
extended pattern of depolarization canals in maps of polarized 
intensity. They constructed an RM map from these five wave- 
lengths and interpret the ring as a purely magnetic flux tube. 



Send offprint requests to: M. lacobelli; iacobelli@strw.leidenuniv.nl 



[ Bernardi et al.| ( |2009| ) observed the FAN region at A = 
1.93 - 2.16 m in order to characterize the properties of the fore- 
ground for epoch of re-ionization experiments. They detected 
total intensity and polarization fluctuations for the first time 
in the Galactic diflTuse foreground emission at this wavelength. 
Rotation measure synthesis (RM-synthesis) was used to image 
this region in the sky revealing a complex and structured distri- 
bution of polarized signals. 
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In this paper, we re-analyse the data from Bemardi et al. to 
study the diffuse Galactic polarized emission. In Section 2 we 
present the data and illustrate the working principle of RM syn- 
thesis technique. In Section 3 we describe the main properties of 
total and polarized intensity. In Section 4 the polarized data in 
Faraday depth space are discussed. In Section 5 we describe the 
main components and their properties. In Section 6 we present 
a simple model to explain the observed structures and their fea- 
tures. Finally in Section 7 we summarize and conclude. 



2. Data analysis 

The data were collected with the east- west Westerbork Synthesis 
Radio Telescope (WSRT) array in 2007. Baselines between 36 m 
and 2.7 km were used to obtain good wv-coverage. Data reduc- 
tion was performed by Gianni Bernardi using the AlPS-h-h pack- 
age, and for a description of their reduction we refer to Bemardi 
et al.| ([2009 ). The corresponding main properties are given in 
Table 1 . Ionospheric propagation effects are direction- and time- 
dependent and affect the phase of the signal received by the inter- 
ferometer. They also affects the polarization, giving rise to time- 
variable Faraday rotation. Moreover, differential Faraday rota- 
tion by the ionosphere will alter Stokes Q and U, which could 
lead to depolarization but could also lead to spurious extra po- 
larization, thus affecting the results of the polarization imaging. 
Over the six days of observation the mean total electron content 
(TEC) value (provided by the Center for Orbit Determination 
in Europe (CODE)) and the mean intensity of the geomagnetic 
field (calculated using the International Geomagnetic Reference 
Field (IGRFIO)) give us a typical RMion < 0.3 rad m"^, which 
implies a phase shift of about 80 degrees at the lower observed 
frequency. The spread of the RMion over days of observation 
is 0.1 rad m"^. Since the resolution in Faraday depth is about 
3 rad m^ (see Table 1), ionospheric Faraday rotation will hardly 
change the analysis or conclusions, so that corrections for the 
ionospheric Faraday rotation were not applied. 
To enhance sensitivity for diffuse emission, the resolution of the 
Stokes Q and U maps was decreased. In this paper we use the full 
angular resolution Stokes I map (i.e. A^o = 2 arcmin) and Stokes 
Q,U maps with angular resolution = 4.2 arcmin. No correc- 
tion was performed for the primary beam attenuation in order 
to preserve a uniform noise level across the images. In the fol- 
lowing section, we summarize the principle of the RM-synthesis 
technique. 



form between F((f)) and the measured complex polarization vec- 
tor P(A^): 



2.1. RM-synthesis 

Low frequency polarimetric studies of the diffuse Galactic radio 
emission are a valuable tool to investigate the ionized and mag- 
netized components of the interstellar medium. Generally, these 
studies have been based on RM measurements calculated as the 
slope of a linear fit to the polarization angle as a function of 
the wavelength squared. However, a new analysis t echnique ha s 
recently been applied, the so-called RM-synthesis ( Burn|1966| ). 
This technique has only recently become practically applicable 
due to technical and computational advances, and was develo ped 
and used for the first time by |Brentjens & de Bruyn|pQ05| ). Its 
main advantages with respect to the standard RM technique are 
that it provides a mapping of the linearly polarized emission as 
a function of Faraday depth (cp), the Faraday dispersion func- 
tion (F((f))), and that bandwidth depolarization and depth depo- 
larization are much less severe. RM-synthesis is implemented 
as a weighted (by the sampling function W(A^)) Fourier trans- 
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where the electron density He and the magnetic field B are given 
in cm"^ and yuG, respectively. Since we can only sample a fi- 
nite positive range of wavelengths, the inversion of the Fourier 
transform is incomplete, and we only get an approximation of 
F(0). However, three physical quantities directly linked to the 
experimental set up (the channel width d??, the width of the A? 
distribution A/1^, and the shortest wavelength squared /l^.^) can 
be used to characterize the detection capability: 

- the maximum detectable Faraday depth, constrained by the 
channel width 6A\ (pmax ~ ^/SA^; 

- the maximum scale detectable in Faraday depth to which the 
sensitivity is reduced to 50%, which is constrained by the 
lowest observed wavelength Amm'- ^(j) scale ~ ^/^L>2' 

- the resolution in Faraday depth space defined as the half 
power width of the rotation measure spread functiorF] 
(RMSF), constrained by the total wavelength range A/1 . 
A0 2 V3/A/l^). This defines the minimum separation be- 
tween separate synchrotron emitting structures that is de- 
tectable. 

For our frequency coverage, ^max ~ 2650 rad m"^, descale ~ 
0.85 rad m"^, and A0 ^ 3 rad m"^. With these data we can there- 
fore only detect Faraday thin structures, since the resolution is 
greater than the maximum detectable scale. Since the polarized 
radiation is a vector quantity, astrophysical information is stored 
in both intensity and angle. However, due to the generally low 
signal-to-noise ratio (S/N) per resolution element, the solution 
of the ^TT-ambiguity is challenging so we do not use polarization 
angle maps for a quantitative derivation, but only to infer some 
basic properties. For this reason, we focus on polarized inten- 
sity maps and maps of Stokes U. Maps of Stokes Q show similar 
characteristics to Stokes U. 



2.2. RM-synthesis data cubes 

Maps for Stokes Q,U, and PI were obtained in Faraday depth 
space by applying the RM-synthesis technique to the 4.2 arcmin 
Stokes Q,U maps in the frequency domain. No bias correction 
was applied to the polarized intensity maps but owi ng to the 
hi gh side lobe level of the RMSF, an RM CLEAN ( ^Heald et 
aL]|2009 ) was performed with a threshold of 5 mJy. The polar- 
ized emission from our Galaxy is ex pected at low Fa raday depth 
values \(p\ < 10^ rad m"^ (see Cleg g et al. 1992 ). Inspection 
of this range of Faraday depth reveals the presence of signifi- 
cant polarized emission only within a few tens of frames around 
= rad m"^ [Bemardi et al.| ( |2QQ9] ). The PI content of the data 
out of this range of Faraday depths has not been checked. For the 
data analysis, we limited the Faraday depth range of this data set 
to [-50, -h50] rad m'^. 



^ The point spread function in Faraday depth space is generally called 
the rotation measure spread function. 
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Table 1. Observational set up and basic properties of the RM-synthesis cube, from Bernardi et al. (2009 ) 



Data 



lb 


137° +7° 


Galactic coordinates 


^bands 


8 


Number of spectral bands 


yO,band 


139.3, 141.5, 143.7, 145.9, 






148.1, 150.3, 152.5, 154.7 MHz 


Central frequency of each band 


^yO,band 


2.5 MHz 


Width of each band 




9.8 kHz 


Frequency resolution 


AOo 


2' X 2.2' 


Beam size 


(a, 6, 0)-cube 


A0 


3 rad m ^ 


RMSF width 


^(l>scale 


0.85 rad m-^ 


Max scale in 


66 


4.2' 


Beam size 




1 mJy bearn"^ = 1 K 


Conversion factor 



3. Observational results 

Two well known basic properties of the diffuse polarized emis- 
sion (see |Wieringa et al.|1993| ) are visible in the images (see e.g. 
Figs.[2]and|8])i 

- the intensity attenuation due to the WSRT primary beam 
(about 2.5° at 150 MHz). This is a strong argument for ex- 
cluding an instrumental origin of the polarization fluctua- 
tions, and it suggests a wider extent on the sky of the detected 
patterns of polarized emission; 

- the absence of a correlation between total intensity and po- 
larized intensity, which indicates that the polarization fluc- 
tuations are caused by Faraday rotation of the difl'use syn- 
chrotron emission background by the intervening magneto- 
ionized medium. 

In the following sections, we describe the global properties of 
the signals detected in polarized intensity and total intensity. 

3.1. Total intensity 

The Stokes I map displayed in Fig. [T] has a size of about 12° in 
declination and 12° in right ascension, and its dynamic range is 
about 1450:1 at an angular resolution of T x 2.2' FWHM. The 
only extended Galactic object detected is the Hii region com- 
plex W3/W4/W5 in the south-west corner, approximately from 
((^,^)/20oo ^ (02^25-, 62°. 0) to (cr,^)y2ooo ^ (02^53^, 60° .4). 
Faint Stokes I emission is also seen from the spiral galaxy IC 342 
at (a,^)/2ooo ~ (03^7"^, 68°. 0) and from the giant double ra- 
dio galaxy WNB0313-h683 at (of,^)/2ooo ~ (03^8'^, 68°. 3). 
No large-scale difliise mission in the primary beam is seen 
at 150 MHz because the interferometer filters away power on 
the largest scales. However, B ernardi et al., (2009) detected ex- 
tended faint, mottled, small-scale total intensity emission from 
the Galactic foreground. From the power spectrum these authors 
suggest a minimum spatial scale of ~12 arcmin for these Stokes I 
fluctuations. 

3.2. Polarized intensity 

Figure |2] shows the polarized intensity at a Faraday depth = 
-1 rad m"^. Grating rings are seen around brighter compact 
sources, but owing to their multiplicative nature, such error pat- 
terns are only visible around bright Stokes I sources, so their 
efl'ect is negligible in the map except in the direct vicinity of 
bright point sources. Unwanted instrumental signals are visi- 
ble in the Stokes Q, U, and PI maps as elongate d stripes and/or 
waves crossing the field. As pointed out by(de Bruyn & Brentjens 




Right Ascension (J2000) 



Fig. 1. Total intensity map at 150 MHz with an angular resolu- 
tion of^Tx 2.2'. The diffuse structures towards the upper left 
and lower right corner are the galaxy IC342 and the Hii region 
complex W3/W4/W5 respectively. 

('2005^ and 'P izzo QtH] (jgOTT), the WSRT off axis polarization 
has a strong frequency dependence with a period of 17 MHz, 
which causes peaks in the Faraday spectrum at about 42 rad m"^. 
At our frequency range, this 17 MHz ripple appears as a broad 
peak at the edges of the selected Faraday depth range (see 
Fig.|4]). These are due to sources Gas A and Cyg A entering the 
telescope side lobes. 

The overall properties of the polarized signal distribution as 
a function of Faraday depth are shown in Figure |4] and an ex- 
ample of Faraday dispersion spectrum is displayed in Figure [3] 
The distribution of polarized signal includes the whole field, 
i.e. 480 X 480 lines of sight. We counted at each Faraday 
depth the number of pixels having a polarized intensity above 
a 5 mJy beam"^ threshold (~ 5 times the PI mean noise level). 
Polarized emission is detected over a wide range from positive 
to negative values. Two main features are clearly seen: 

- a central peak at slightly negative cp where the most of polar- 
ized emission is. 
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Fig. 2. Frame of the RM-synthesis cube at Faraday depth = 
- 1 rad displaying the extended polarized emission. Grating 
rings can be seen around the bright point sources at {a, 6) /2000 ~ 
(02^38'^, 59° . 1) f close to the edge of the image), and {a, 5) /2000 ~ 
(03^06^, 62°. 5). The drawn boxes are chosen to infer the noise 
behaviour (box 1) and derive the cross -correlation profile (box 
2 ), see text. 
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Fig. 3. Example of Faraday dispersion spectrum with emission 
at a single Faraday depth along the line of sight, averaged over 
a beam-sized box region. The noise level is indicated by Icr error 
bars. 



- two wings towards the edges of our range, due to the resid- 
uals from Cas A and Cyg A. 

Figure |4] makes it evident that the Galactic emission from this 
region is concentrated in a very narrow (p range from about -10 
to +5 rad m"^. The main central feature in the distribution is 
clearly asymmetric; we interpret this as an evidence of a multi- 
component signal and we discuss this in Sects. 5 and 7. 



-40 -20 



20 40 60 



;> [rad m~-] 



Fig. 4. Distribution of polarized intensity as a function of 
Faraday depth, above a threshold of 5 mJy. At |0| > 40 rad m"^ 
the residual signal from side lobes of Cyg A and Cas A are visi- 
ble, while the Galactic emission occurs around -1 rad m~^ and 
has an asymmetric profile. 

4. The polarization cube in Faraday depth space 

In this section we discuss the noise properties and the connection 
between RMSF width and cross-correlation in Faraday space. 



4.1. Noise properties and errors 

Since we are dealing with the measurement of a pseudo- vector - 
i.e. the polarization - in the presence of random noise we look at 
noise behaviour. While the two Stokes Q and U parameters are 
Gaussian distributed, the linearly polarized intensity 



PI= Vg^ + 



(3) 



follows an asymmetric positive definite distribution: the Rice 
distribution, which is dependent on signal-to-noise and respon- 
sible for a bias towards too high values of PI at low S/N Wardle 
|& Kronberg| ( [T974l ). 

We tested the Gaussian behaviour of Stokes Q and U maps 
and estimated the noise in PI as a function of Faraday depth un- 
der the assumption of a Rayleigh distribution of polarized inten- 
sities. The noise was determined from a small region assumed to 
contain no emission from the sky, so that all detected emission is 
noise (see box 1 in figure [2]). We assume uniformity of the noise 
across the field. Since in every Faraday depth frame we detect 
emission (either from the ISM or from grating lobes from strong 
out-of-field sources), this uniformity could not be tested. Noise 
values in Stokes Q and U are similar; also, the assumption of 
Gaussian noise in Stokes components Q and U is satisfied, be- 
cause the standard deviation of the sample and the width of fitted 
normal distribution agree with each other. 

The noise properties of PI as a function of the Faraday depth 
under the assumption of Rayleigh distribution statistics are 
shown in figure [5] Good agreement is found for the expected 
(i.e. cTpi = 0.66(T[/) and the measured noise variance, as well as 
for the expected (i.e. < /dpi >= 1.25(T[/) and the measured noise 
mean value. A mean noise level of ^ 1 mJy beam"^ is found 
in PI over a range of [-50, -i-25] rad m"^ while an increase 
in the mean and variance values is seen in the Faraday depth 
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Fig. 5. Noise behaviour in the PI map derived from a region 
denoted by box 1 in figure^as a function of Faraday depth un- 
der the assumption of Rayleigh distribution statistics. The ex- 
pected mean and variance values of the noise are indicated by 
the dashed and solid lines, while the measured values indicated 
by circles and stars refer to the noise mean and variance, re- 
spectively. 



range above +25 rad m"^ due to the presence of residual signal 
from side lobes of Cyg A and Cas A. Because the Rayleigh 
statistics a clipping of PI values at Scr allows to mitigate the 
impact of Ricean bias below 10% level. Since our analysis does 
not depend on exact values of polarized intensity, but rather uses 
rough estimates of ISM parameters, bias corrections were not 
applied. 

No corrections for the off-axis instrumental polarization 
were applied. We estimated its magnitude and radial dependence 
from the fractional polarization values at = rad m"^, where 
instrumental polarization is expected. We found a mean level of 
^ 2% within 3° (primary beam), which rises to ^ 9% at the 
edges of the field of view. There are not enough sources left in 
the maps after the subtraction to derive the angular dependence 
of the antenna pattern. 

In analogy to radio interferometric observations, the standard 
errors in Faraday depth measurements can be obtained in an RM- 
synthesis cube as: cr^ = 0.5A(f)/(S /N) ~ 1.5/(^/A^) rad m'^ 
where A0 is the RMSF width and S/N the PI signal-to-noise ratio. 
Because we mostly deal with spatially extended emission over 
several ranges of Faraday depths, we prefer to provide a general 
upper limit for the cp uncertainty at each frame as set by the mean 
S/N of the emission across the maps. Because the chosen 3cr clip 
level (i.e. S /N > 3.0): cr^ < 0.50 rad m'^ Therefore the total 
Faraday depth estimation error in our 2 m data consists of two 
main contributions: the systematic error due to the ionospheric 
Faraday rotation (as discussed in Section [2]) and the standard 
error (as shown in Figs.[3]and[9]), having comparable amplitudes. 

4.2. Cross-correlation in cp space 

In RM- synthesis the standard estimator for the nominal reso- 
lution in Faraday space is the RMSF width. However, as in 
aperture synthesis, there is some information present on scales 
smaller than the beam size, as can be seen in figure |7] This 
figure shows polarized intensity changes between neighbouring 



frames (i.e. A0 = 1 rad m"^), about one third of the RMSF 
width. An eff'ective approach to investigating 0-resolution is 
cross-correlation. Cross-correlation between frames can be used 
as a probe to infer the presence of structures in the data and their 
Faraday thickness. We calculate the cross-correlation coefficient 
for the frames a and b with a frame gap Sep (Cs(p) as 



1 ^ (ciij - a)(bij - b) 



^ 

O-aO-b 



(4) 



where N is the total number of pixels, (a, b) are the mean val- 
ues of the frame, and (cr^, crij) are their standard deviations. The 
frame gap Sep = \a-b\i^ the distance between frames a and b in 
rad m"^. 

We studied the correlation length in Faraday depth space of 
the PI data by considering the whole field of view and by fo- 
cussing on a square region in the centre of the field (given by 
box 2 in figure [2]). In figure [6] we show the cross-correlation co- 
efficients for diff'erent frame gaps (6^). 

For a fixed convolving RMSF, the reference level for data 
correlation/anti-correlation is shifted from zero to a certain pos- 
itive value. This off'set is due to the convolution of the noise of 
each frame with the RMSF, which introduces a degree of corre- 
lation of the noise in the frames. The cross-correlation of frame 
gaps d(p = 1, 2, 3 rad m"^ for noise-dominated frames convolved 
with an RMSF of 3 rad m"^ is approximately 0.68, 0.22 and 0.1, 
respectively. Therefore the average levels in Fig. [6] correspond 
to a convolving function with a width of ^ 2.95 rad m"^, in 
good agreement with the width of the RMSF. This means that 
outside the range of Galactic signal of the Faraday depth, from 
about -10 to +5 rad m"^, the cross-correlation profile is con- 
sistent with a pure noise signal. Therefore for S(p = I rad m"^, a 
correlation coefficient above 0.68 indicates a positive correlation 
between two frames, a correlation coefficient of about 0.68 indi- 
cates no correlation while a correlation coefficient below 0.68 
indicates an anti-correlation. 

Several features in the cross-correlation profile for Sep = 
1 can be seen, which disappear as the frame gap increases. 
Therefore they are not broader than the RMSF width, which 
means these structures are Faraday thin. For a frame gap of 
Sep = I rad m"^, prominent peaks are seen at -5 rad m"^ 
and (p ^ -2 rad m"^. Three more tentative, small peaks, pos- 
sibly related to fainter extended emission, are observed around 
(p ^ -10 rad m"^, 0^+1 rad m"^, 0^+5 rad m"^. 

Calculating the cross-correlation coefficients with the same 
frame gap of ^0 = 1 rad m"^ in the smaller box indicated 
as 2 in Figure |2] shows interesting results. Correlated peaks at 
-10, -6, -2 rad m"^ are still seen, although their amplitudes 
have changed. However, the most conspicuous changes are seen 
in the range of Faraday depth ~ [-2, +4] rad m"^ where the 
cross-correlation profile is turned into anti-correlation, with two 
possible small peaks at rad m"^ and (p ^ -\-2 rad m"^. The 
polarized emission in this 0-range is spatially anti-correlated. 
Simple simulations show that qualitatively this behaviour can be 
caused by a spatial gradient in Faraday depth space across these 
frames. The presence of multiple Faraday components within the 
peak was also confirmed by comparing the width of the peaks in 
the spectrum with the width of the RMSF. We checked the be- 
haviour of some representative lines of sight through the "bub- 
ble", "ring", and "curtain" components. From the fitting of the 
limited sample of Faraday dispersion spectra, we find clues of 
multiple poorly separated components, in agreement with the 
cross correlation result. 
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Faraday depth [rad m -] 

Fig. 6. Cross-correlation coefficients as function of Faraday 
depth for different frame gaps Sep = I rad (solid), Sep = 
2 radm~^ (dash-dotted), and Sep = 3 radm~^ (dotted). The cross- 
correlation at S(p = 1 rad for only the emission within box 
2 in figure |2] is shown by the dashed line. The uncertainties are 
shown by \cr error bars. 



- the RM peaks associated with the noise or artefacts (e.g. 
due to the lobes of the RMSF) show a systematically low 
FWHM < 2.7 rad m"^, while RM peaks associated with 
the three components have FWHM ~ [2.9, 3. 2] rad m"^, in 
agreement with the expected RMSF width; 

- Some lines of sight show excess of widths or complex mor- 
phology (see e.g. spectrum of line of sight 4 in Figure [9]), 
which make them unsuitable for a single component fitting 
step. These cases can be found in all three components, but 
are more evident for the "curtain" component; 

- Often spectra also have associated main lobes that are higher 
than the ~35% level expected from the RMSF profile. In con- 
clusion we find agreement between the result of this test and 
the cross-correlation. 



5. Definition and description of structures 

The (a, S, 0)-cube contains a wealth of structures. In this section 
we briefly describe the main morphological features in the PI 
and Stokes U maps as a function of Faraday depth. First, there 
are spatially compact and isolated structures, such as Galactic 
and extragalactic Stokes I sources. These objects primarily "il- 
luminate" the intervening Galactic ISM, and reveal its Faraday 
depth. If their structure is multi-component, then internal 
Faraday dispersion can be significant and they may actually 
show emission over a range of Faraday depths. Their definition 
is based on their compactness in (a,S) and, occasionally, on 
association with known objects. The second class of structure is 
the spatially (very) extended emission detected in polarization. 
These can be characterized by their properties of PI, Q, or 
U in Faraday space, in combination with the distribution of 
polarization angles. 

We identified various components in the RM-synthesis cube 
based on morphological consistency and coherence in Faraday 
depth and name these components the "bubble", the "ring", and 
the "curtain". These components are shown in polarized inten- 
sity in Figure [8] 



- The discrete PI emission from the "bubble" appears at ~ 
-10 rad m"^ in the centre of the field of view. A circular 
filled structure develops and expands with a maximum of 
polarized intensity at 5 rad m"^ and an angular diam- 
eter of -3°. Also, close to the center a hole in the emission 
develops and expands until the next structure - the "ring" - 
comes up. The Stokes U map (see Fig. [7] upper left panel) 
depicts a complex pattern of small-scale wrapped structures 
surrounded by larger and smoother ones. 

- The "ring" structure consists of extended emission combined 
with narrow canals of low polarized intensity around a cen- 
tral, almost circular, edge. The "ring" can be considered as a 
smooth transition from the "bubble". Its polarized emission 

reaches a maximum around (p 1 rad m"^ and shows bright 

PI emission occurring in the north-western part of the field 
with longer depolarization canals ([ Haverkorn et al. 2000), 
with lengths up to a few degrees. The orientation of these 
canals is mostly paralle l to the Galactic plane, i n agreement 
with the observation of |Hayerkom et al.| ( |2003| ). The Stokes 
(Q and) U maps (see Fig. [7j upper right panel) show a rapid 
sinusoidal behaviour at the edge of the "ring", corresponding 
to a gradient in polarization angle. 

- Finally, a significant pattern of faint, extended polarized 
emission -the "curtain"- is found at ~ [0, -h5] rad m"^. A 
fast transition is seen around (p ^ -\-2 rad m"^ from the pre- 
vious "ring" component to the "curtain". We interpret this 
as a spatial gradient of polarized emission in Faraday space 
around ~ -2 rad m"^ and ~ rad m"^, which is re- 
sponsible for the abrupt decrease in correlation. Stokes U 
emission (Fig. |7] lower panels) reveals extended features on 
larger scales than the "bubble" and "ring" components. Also, 
whereas the "bubble" and "ring" show morphological corre- 
spondences, these are not present for the "curtain". 



These three components in the RM-synthesis cube correspond 
to the peaked features seen in the cross-correlations profiles 
(Figure]6]) around 6, -2, -\-2 rad m"^, respectively. 



5.1. Diffuse foreground polarized emission 

To emphasize properties of linearly polarized emission in 
space, the PI emission was integrated in the frames from 
-13 rad m"^ to -5 rad m"^, from - 4 rad m~^ to - 1 rad m"^, 
and from rad m"^ to -h5 rad m"^ as (Brentjens 2007 ) 



(pmax j — 



(5) 



where the normalization factor B is the area of the RMSF beam 
after RM-CLEAN divided by the interval between two frames 
of the data cube, and crpj is the noise in PI estimated from 
figure [5] These maps were clipped at a level of 3cr and then 
combined into a composite false-colour image as shown in the 
left-hand upper plot in figure [8] 

Figure |9] shows a sample of Faraday depth spectra, obtained 
averaging over a beam-sized box region, for several lines of 
sight. The lines of sight were selected to show widely varying 
spectra, with a strong main peak (e.g. spectrum 7) or with dou- 
ble and multiple peaks (e.g. spectra 1, 3, and 5 respectively). A 
peak in the Faraday depth spectrum indicates an emission com- 
ponent in Faraday space. However, due to the limited frequency 
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ghggm 3^00^" 2^30^^ 3^30"" 3^00"" 2^30"^ 

Right Ascension (J2000) Right Ascension (J2000) 

Fig. 7. Stokes U emission at four representative Faraday depth values. The images correspond to (p = -5 rad m~^ ( upper left panel), 
(p = -1 rad m~^ (upper right panel), = rad m~^ (lower left panel) and = +1] rad m~^ (lower right panel). 

Table 2. Properties of the three components. 



Structure 


Colour 


range 
[rad m~^] 


4^ peak 

[rad m~^] 


<P> 
[Jy bearn"^] 


Bubble 


Red 


8.0 


-5 


(1.72 + 0.58) X 10"^ 


Ring 


Green 


4.0 


-2 


(1.78 + 0.26) X 10-3 


Curtain 


Blue 


5.0 


+2 


(1.53 + 0.22) X 10-3 



Notes. The columns denote: 1) Name of component, 2) Colour in figure |8] 3) Range in Faraday depth in which component is visible, 4) Faraday 
depth of peak, 5) Polarized surface brightness averaged over Faraday depth range. 

coverage, components in Faraday space may be represented in higher than the "curtain". The mean polarized brightness tem- 
the Faraday depth spectrum incompletely or not at all. perature over the range g [-13, +5] is ^ 5.6 K. 

We estimated the mean polarized surface brightness for each 
component (see Table 2). These mean values are of the same 
order of magnitude. Interestingly the "bubble" and "ring" com- 
ponents have comparable polarized surface brightness, slightly 
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Right Ascension (J2000) 




- o.ooao ^ 
33 



Right. Ascension 




O.OOiJO ^ o 
2q 



3''3Q™ 3^00^ 
Right. Ascension 




3''3Q™ 3^00^ 
Right. Ascension 



Fig. 8. Composite image showing in the upper left panel a colour coded image of polarized intensity emission clipped at 3cr at 
three main Faraday depth ranges, where the white crosses denote the lines of sight used to extract the Faraday dispersion spectra in 
figure^ PI clipped maps averaged over ranges in Faraday depth depicting the ''bubble", the ''ring" and the "curtain" are shown 
clockwise in grey scale; Red ( upper right panel) is averaged PI over a range ofcp e [-13, -5] rad m~^, Green (lower right panel) is 
averaged PI over a range ofcp e [-4, -1] rad m~^ and Blue (lower left panel) is averaged PI over a range ofcp e [0, +5] rad m~^. 



5.2. Polarized extragalactic background sources 

Although we are primarily interested in the Galactic emission we 
briefly point out a few interesting discrete extragalactic objects 
within 5° of the field centre. These are 

- The spiral galaxy IC 342, which is visible both in polarized 
and total intensity. Large scale disk/halo emission at (/) = 
(-3.0 ±0.5) rad m"^ is observed, in agreement with Krause 
et al. |l989| ) who find RM values between -10 rad m"^ and 



10 rad m . 

- The giant double-lobe radio galaxy WNB 0313+683, which 
shows prominent linearly polarized emission from the 
south lobe (at (a,^)/2ooo ~ (03^7^^, 68°. 2)) and the 
core (at (a,S)j2ooo ~ (03^8^^, 68°. 3)) in a range ^ 



ISchoenmakers et al.| ( [1998] ) and [Haverkom otsl] ( [2003] ). 
However, no polarized emission is detected from the north 
lobe, which is pointing away from the observer. Since the 
Stokes I map shows a sharp diff'erence in brightness for the 
lobes (the north to south lobe mean brightness ratio is -0.2), 
the non-detection of polarized emission from the north lobe 
can be explained as a consequence of relativistic beaming 
and/or additional Faraday depolarization (Heald et al. 2009 ). 

We used discrete polarized background sources to search for a 
latitude and/or longitude dependence of the Faraday depth across 
the field of view. From the NVSS RM catalogue ( Taylor et al. 
2009 ) we selected a sample of polarized sources with known 
RM, which covers the observed field of view and its surround- 



[-17,-5] rad m ^ consistent with values reported by ings (see Fig. 10). To average out intrinsic RM components of 
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Line of sight 1 Line of sight 2 Line of sight 3 





Fig. 9. Faraday dispersion spectra of selected lines of sight labelled from 1 to 9 as in figure^ The three dashed black lines point 
out the position of the histogram peak for the ''bubble", the ''ring" and the "curtain" components. For each spectrum the "Clean" 
line is obtained from the "Raw" one by applying an RMCLE AN session with a threshold of 3 mJy beam~^. The noise level is shown 
by \cr error bars. 



the extragalactic sources, we re-sampled the data over square 
bins with an angular scale of 5° (about half the size of the 
imaged field). We find a clear gradient in Galactic latitude of 
ARM/Ab ^ 4.6 rad m"^ deg"\ while no evident Galactic longi- 
tude dependence is found. The dominant negative RM values 
imply overall negative B\\ values, in accordance with the ob- 
served negative Faraday depths. A systematic increase in RM 
values for sources below b = 5° is observed as expected be- 
cause of the presence of the Galactic disk. In the longitude range 
/ G [125°, 140°] the presence of extended Hii regions (i.e. the 
W3/W4/W 5/HB3 H ii region/supernova remnant (SNR) complex 
in the Perseus arm), as well as the Ha maps ( IHaffner et al.|2Q03| 
suggests that the increase in RM towards lower latitudes is at 
least partially a consequence of enhanced thermal electron den- 
sity. We note that both the diff'use emission and polarized point 
sources show predominantly negative Faraday depth and RM, 
respectively. However, RMs from the extragalactic point sources 
are much more negative than the Faraday depths of the diff'use 
emission, suggesting that (a) the diff'use emission is nearby and 
does not span the full line of sight through the Galaxy, and (b) 
that |RM| increases through the outer Galaxy without large-scale 
magnetic field reversals. 



6. Model for the diffuse polarized emission 

In this section, we discuss plausible sources in the ISM that 
correspond to the three polarized components identified in 
Section 5. Interpretation of these polarized components in terms 
of synchrotron emitting and/or Faraday rotating regions is not 
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Fig. 10. Distribution of polarized extragalactic point sources 
and rotation measure values from the NVSS catalogue after bin- 
ning over 5 degrees spatial scale across the imaged field and its 
surroundings. The position of pulsar B0226-\-70, which is used in 
section 7 to constrain the spatial distance of observed polarized 
structures is indicated by a cross. 



straightforward for a number of reasons. First, these are inter- 
ferometric data, which are missing short spacings in the uv- 
plane and are therefore insensitive to structure on scales larger 
than two degrees. This explains why a large amount of polar- 
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ized emission is obs erved, but only very weak total intensity 
( Bernardi et al.|2009] ), providing an apparent polarization degree 
far above 100%. This indicates that the synchrotron radiation is 
emitted on large scales, but its polarization is altered on smaller 
scales by Faraday effects. In this case we can calculate the miss- 
ing large-scale structure contribution for Stokes Q and U for a 
uniformly polarized background propagating through a small- 
scale Faraday screen ( Hav erkorn 2 002 ). The expected offset for 
Stokes Q and U (Qq and Uq) depends on the wavelength (A) and 
the width of the Faraday depth distribution (cr^) as 







Qo = Uo = 0. 



(6) 



At low frequency this condition is typically well satisfied and 
we do not expect large-scale undetected Stokes Q and U struc- 
tures. Although there is missing large-scale structure in total 
intensity Stokes I, there is no missing large-scale structure in 
Stokes Q and U, therefore the polarization data are reliable. 
Second, the limited frequency range of the data set precludes 
the detection of Faraday d epth components that have a width 
(50 > 1 rad m"^ (Section \2A\ . Therefore, by definition, the 
Faraday depth components that we observe are unresolved in 
Faraday depth. Analogously to the missing short-spacing prob- 
lem, this indicates that there likely is Faraday-thick polarized 
emission (e.g. emission over a wide range of Faraday depths), 
which we cannot detect. Last, depolarization effects can decrease 
the polarized intensity of the observed structures. 

First, we calculate some properties of the global interstellar 
medium in the direction of the FAN region, below. Then we dis- 
cuss what information is available on the absolute and relative 
distances of the three components in Section [6?T] Finally 



present our model in Section 6.5 



we 



From the mean dust-corrected Ha surface brightness across 
the (half upper part of the) field of view (I Ha ~ 3.3 R) and 
adopting an electron temperature in the diffuse ionized interstel- 
lar medium T = 8000 K, we derive the related emission measure 
EM field = 7.43 pc cm"^. This corresponds to a large-scale mean 
free electron density of 3 x 10"^ cm"^ over an assumed 
total path length o f 7 kpc, a value, which is in good agreement 
with the literature ( |Hill et al.|2008l |Ferriere|200T| |de Avillez et| 
|al.|2012| ). 

The Fan region has a high degree of polarization, which in- 
dicates a relatively regular, large-scale magnetic field topology. 
This can be either uniform or anisotropically random magnetic 
fields. 

The uniform magnetic field component following the spiral 
arms is almost perpendicular to the line of sight in this direc- 
tion. Therefore, the impact of the large-scale interstellar mag- 
netic field on Faraday depth values toward the Fan region is ex- 
pected to be small. We can test this assumption by estimating the 
path-averaged magnetic field amplitude along the line of sight up 
to a few hundred parsecs. We consider a sample of nearby pul- 
sars around the imaged field (see Table [3]) and apply the general 
relation: 



^ DM 



(7) 



With the exception of object 1, all pulsars have estimated 
distances lower than ^ 600 pc. The scatter in B\\ is high, and its 
sign changes for different pulsars, confirming that the uniform 
magnetic field component is much smaller than the random com- 
ponent in this direction. We calculate a mean parallel magnetic 
field {B\\} ^ 1.2 yuG. 



6.1. Distance estimates to the structures 

Most of the polarized emission can be found at slightly negative 
Faraday depths, consistent w ith earlier observations (Brou w &| 
|Spoelstra|1976 [Taylor et al.|2009 ), which suggests that the radi- 
ation is emitted at a short distance from th e observer. IWilkinsonl 
& Smith] (| 19741) and'Uyaniker et al.|(|2003 ) deduce that this po- 



larized emission and the low RM values in the second quadrant 
around / ~ 150° must be nearby, possibly at a distance < 500 pc. 
Also, the NVSS gradient show a mean RM value at the FAN 
latitude of RMiy^qo ^ -32 rad m"^, which is higher than the ob- 
served Faraday depths. We conclude that the detected structures 
are not representative of the general behaviour through the entire 
line of sight; as a consequence, these structures must be close. 

We can roughly estimate the distance to the "bubble", as- 
suming uniform conditions for the medium along the entire line 
of sight, as 



, (pOBS ^ r 

do = D RE pipe] , 

(pREF 



(8) 



where (poBS is the measured Faraday depth, and (Pref the 
reference value of Faraday depth (or rotation measure) cor- 
responding to the line of sight total length Dref- We con- 
sider the rotation measure of nearby radio pulsar B0226-h70 
at (q^, ^)/2ooo ~ (02^31"^, 70°. 26) as a reference. The distance 
of this pulsar is dpsR = (2.25 ± 0.56) kpcn and its rotation 
measure RM psr = (-56 + 21) rad m ~^ from the ATNF Pulsar 
Databas^ ( |Manchester et al. 2005 ). Equation (8) then gives 
dBUBBLE ^ (201 + 90) pc, whic h is consistent with a previous 
estimate ( Haverkorn et al.|2()Q3] ). At such a distance and consid- 
ering the angular size of the "bubble" 3°), a physical size in 
the plane of the sky of 10 pc is obtained. 

Considering the intrinsic Faraday depth for the "curtain" and 
"ring" components along with the large-scale electron density 
field and magnetic field component estimated above along the 
line of sight, we obtain an estimate of the distance needed to 
build up this Faraday depth. We find ranges of ^ 70 pc and 
^ 140 pc for the "curtain" and the "ring", respectively. These 
distance estimates do have considerable error margins due to the 
strong assumptions. 

These distance estimates can be tested by comparing the 
amount of polarized intensity built up over that distance to the 
observed polarized intensity in the component. At 408 MHz, 
the total inten sity emissivity is e^^^ - (11 + 3) K kp c~^ or 
0.01 1 K pc"^ ( jBeuermann, Kanbach & Berkhuijsen|l985| ). With 
a spectral index of -2.5, the conversion factor from 408 to 
150 MHz is ~ 12.2. Thus the total intensity emissivity at 

(134 ± 37) K kpc-i or 



150 MHz is 
0.13 K pc- ^ 



U50 



12.2 X e\ 



.408 



(|Haverkorn et al 
~ 0.0094 K pc 



8% 

2008k, then ei]^ ~ e]^^ x 0.08 = (10.72 ± 



If we as sume a polarization degree of 



-PI 

^ Therefore, 



for a constant 6^^^ 



2.96) K kpc 

(which implies a constant polarization degree of 7%) over the 
line of sight through the "curtain" component, we need about 
(140 ± 39) pc to build up the polarized intensity of the "cur- 
tain" of ~ 1.2 K. This estimate is higher than our estimate 
above, probably because the "curtain" component is likely to 
show enhanced synchrotron emission in a narrow emission re- 
gion (see Section [6^ . The "ring" component would need about 
(168+46) pc to build up the proper polarized intensity, consistent 
with the estimate above. 

We can also give information about the relative positions and 
distances of the components. The "curtain" component is the 



http ://w w w. atnf .csiro . au/research/pulsar/psrcat/ 
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Table 3. Parameters of nearby pulsars located in the direction of the FAN region. 



Object 


Name 


l,b 


DM 


RM 


Dist^ 


Dist DM^ 


(B\\) 






[deg,deg] 


[cm~^ pc] 


[rad m~^] 


[kpc] 


[kpc] 




1 


B0450+55 


152.62,7.55 


14.495 + 7 


10.0 + 3 


1.19 


0.79 


+0.8555 


2 


B0655+64 


151.55,25.24 


8.771+5 


-7.0 + 6 


0.48 


0.48 


-0.9896 


3 


B0809+74 


140.00,31.62 


6.116 + 18 


-11.7 + 13 


0.43 


0.33 


-2.3721 


4 


B1112+50 


154.41,60.36 


9.195 + 8 


3.2 + 5 


0.54 


0.54 


+0.4315 



Notes. Fo r each object the inferre d magnetic field component along the line of sight (B\\) is also listed. Data are taken from the ATNF Pulsar 

Database ([Manchester et al.|2005] ). 

Dist is the best estimate to the pulsar distance (kpc). DistDM is the distance based on the Taylor & Cordes (|1993|) electron density model. 



nearest component. Since it does not show any imprint of the 
"bubble", it has to be located in front of it. The observation of 
small Faraday depths in this component suggests that it is the 
nearest layer of ISIM with respect to the observer. 

The Faraday depths corresponding to these three compo- 
nents are in adjacent ranges, suggesting that their relative dis- 
tances are small. In particular, the "bubble" has to be located 
close behind the "curtain" . This is because an additional syn- 
chrotron emission component would build up in a large distance 
between "bubble" and "curtain", which would have been observ- 
able. The morphology of PI patterns, as well as the very similar 
mean polarized surface brightness (see Table 2), point to a clear 
link between the "bubble" and the "ring". 

6.2. The foreground "curtain" component 

The change in sign of Faraday depth between the "curtain" 
and the other components indicates a reversal of the magnetic 
field component along the line of sight B\\ behind the "cur- 
tain". A plausible magnetic field configuration of how this can 
be achieved is, e.g., a localized, ordered, magnetic field in the 
Galactic plane perpendicular to the observer, which fans out lo- 
cally to produce oppositely directed parallel field components. 
This magnetic field configuration has a maximum synchrotron 
emission at the location of the perpendicular magnetic field lines, 
possibly corresponding to the synchrotron emission from the 
"curtain". 

We propose that this component corresponds to the wall 
of the Local Bubble. The above size estimates for this fore- 
ground component roughly agree with the dimension of the 
Local Bubble in this direction ( Sno wden et al.|19 98 ). A tomog- 
ra phy of the interstella r gas within 250 pc was recently obtained 
by Welsh et aLl ( 2010| combining interstellar absorption data for 
stars with a Hipparcos measured parallax. Maps of spatial dis- 
tribution for the neutral gas and neutral plus ionized gas compo- 
nents were derived from Nai, Can tracers, displaying the Local 
Bubble as a region that is poor in neutral gas (and dust) but with 
several diffuse and highly ionized clouds. A visual inspection of 
both Nai and Can galactic projected maps (see their Figs. 12 
and 15) toward the direction / = 137°, b = -h7° reveals the pres- 
ence of an extended and high-density wall of partially ionized 
medium at about 100 pc of the Sun, as well as a rarefied medium 
inside the wall, in agreement with our estimate. IMoreover, opti- 
cal starlight measurements in the Fan region from Heiles ( 2000 ) 
show a rise in polarization degree from a distance ~ 100 pc, as 
shown in Fig. 1 1 This can be explained by an increase in the dust 
density, denoting the edge of the Local Bubble. The few high po- 
larization degrees at ~ 10 pc can be explained by a compact and 
cold cloud in the very local interstellar medium. In this picture, 
the magnetic field would be aligned with the Local Bubble wall, 
i.e. mostly perpendicular to the line of sight, which indicates en- 
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Fig. 11. Increase in the polarization degree of optical starlight 
from stars within the Fan field, closer than 3 kpc. 



hanced synchrotron emission in the Bubble wall. The Faraday 
rotation of that emission by the Local Bubble interior is small, 
i.e. < 2 rad m"^. Using the computed mean parallel magnetic 
field component of {B\\) ^1.2 [iG and a path length of 70 pc, 
the mean electron density is n^^ 1.5 x 10"^ cm"-^ in the Local 
Bubble, which agrees with literature ( |Spangler|2Q09| ). 

Previously, an association of the extended polarized Galactic 
ra dio synchrotron fo reground with the Local Bubble was made 
byl Brentjens ( 2011| ), who observed the area around the Perseus 
cluster near / = 150°, b = -13° at .1 = 0.80 - 0.88 m. In this 
region, negative Faraday depth values are expected from the in- 
ferred direction of the large-scale magnetic field, while he finds 
a positive Faraday depth component at = -h6 rad m"^ and esti- 
mates its size of approximately 200 pc. These values diflTer from 
our estimations by about a factor two and may indicate density 
variations and irregular shape of the Local Bubble wall. 



6.3. The background "ring" component 

The "ring" component shows a small south-east to north-west 
gradient in Faraday depth (< 1 rad m"^ degree"^). This gradi- 
ent does not correspond to t he large-scale grad ient in RM ex- 
pected from the NVSS data ( [Taylor et al.|2QQ9| ), but is directed 
almost perpendicular to that. If, as we argue below, the "bubble" 
component is just an additional Faraday rotation of synchrotron 
emission in the "ring" component, we would expect this gradi- 
ent across the field to also be present in the "bubble". However, 
since the Faraday depth gradient of ~ 2 rad m"^ is measured over 
the entire imaged field, the expected gradient over the bub- 



11 



lacobelli et aL: Rotation Measure Synthesis at 2 m wavelength of the Fan region: Unveihng screens and bubbles. 



ble size (about 1/3 of the entire field) turns out to be less than 
1 rad m"^ degree"^ (i.e. less than 1/3 of RMSF width), hence 
undetectable with our observational set up. 

6.4. The "bubble" component 

The "bubble" component is detected in polarized intensity, but 
not in total intensity. This could be due to the lack of short 
spacings. However, if the bubble were significantly emitting, its 
polarized intensity would be higher than the ring component. 
Instead, their polarized brightness temperatures are compara- 
ble, suggesting that the "bubble" is a Faraday rotating structure, 
retro-illuminated by the synchrotron emission from the "ring". 

We interpret these observations as a consequence of the con- 
figuration with the "bubble" located in the spatially extended 
"ring" structure. Two configurations are possible for the "bub- 
ble": it can be assumed to be either a filled sphere or a shell 
surrounding a density depletion. We prefer the first configuration 
because PI emission at the greatest Faraday depths is observed in 
the centre of the "bubble". Indeed under the assumption of pres- 
sure balance a depletion of warm ionized medium relative to its 
surroundings may imply an enhancement of neutral gas or dust, 
which is not observed (e.g. the total H i or IRAS 100 fim maps). 
Instead since the early work of Verschuur] (| 1 969|), a considerable 
lack of H I is observed, as well a very good agreement between 
higher RM / Faraday depth values and lack of H i. Also, for a 
filled sphere it is easier to build up the higher Faraday depths 
observed through the "bubble" 

The additional Faraday rotation in the "bubble" component 
can be caused by increased B\\, increased or a combination of 
these two options. Observations of Ha emission measures set an 
upper limit for the increase in rie in the bubble. The bubble is not 
detected (Haver kom et al.||2003| ) in the WHAM Northern Sky 
Survey (Ha ff^ner et al.|2003| ), which has a sensitivity of 0.05 R. 
For the emission not to be detected at 3cr (i.e. 0.15 R), the emis- 
sion measure has to be EMi,cr ^ 0.34 pc cm"^. 

Assuming that the bubble is spherical, the magnetic field 
component along the line of sight can be estimated as 



^(pmax,Bubble ~ 0.81 X He^Bubble X B\\,Bubble X ^^-Bubble • 



(9) 



We can express both the size and the free electron density as 
a function of the distanc e to the "bubb le": ALBubbie - tan 3° x 
dsubbie and Ue^Bubbie = ^JEM I Msubbie- By using these expres- 
sions we can check the variation in the magnetic field compo- 
nent along the line of sight, the free electron density and the 
size of the "bubble" as a function of its distance, indeed indi- 
cating an enhancement (with respect to the mean values com- 
puted above) of both the free electron density and magnetic field 
strength along the line of sight by a factor of a few as shown in 
figure [12] A reference distance of about 200 pc implies a path 
length through the "bubble" of < 1 1 pc and the corresponding 
electron density and magnetic field strength within the bubble 
feature are Ue^Bubbie -0.18 cm"^ and B\iBubbie ~ 5.25 yuG, respec- 
tively. 

Extra morphology information is obtained from the polariza- 
tion angle maps. A highly structured configuration for the pro- 
jected field is observed in the g [-10, 0] rad m"^ range. As can 
be seen from figure 13 the polarization angles (also in Stokes 
Q,U) show very ordered structures around the edge of the bub- 
ble. The circular, rapid variations in angle are due to a large angle 
gradient, displayed in a fixed angle range of [-0.5;:, 0.5;:] rad. 
Each "ripple" corresponds to an angle gradient over ~ 8 arcmin, 
followed by an abrupt jump across ~ 2 arcmin, which is due to 




160 180 200 
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Fig. 12. Variation in B\\^Bubbie^ ^eMbbie^ (^^d ALBubbie a function 
of ''bubble " distance. Also, for a reference distance of ~ 200 pc 
(vertical dot dashed black lines) the upper limit for the ''bub- 
ble" electron density (ne^Bubbie ~ 0.18 cm~^), the magnetic field 
component along the line of sight B^Bubbie ~ 5.25 jiG and its size 
ALBubbie ~ 10-5 pc are indicated by the horizontal black dashed 
lines. 



the ^TT-ambiguity. The polarization angle thus varies by n radi- 
ans over ~ 10 arcmin. On average four or five of these ripples 
are visible, which indicates an angle gradient of ~ 5;: rad 
over 50 arcmin. At a wavelength of two metres, this indicates a 
gradient in Faraday depth A0 = AO/A^ ^3.9 rad m"^ over 50 ar- 
cmin, corresponding to 2 pc at the estimated distance of the 
"bubble". 

Moreover, we can use the Faraday depth gradient derived 
from the rippled features in the polarization angle map to esti- 
mate the maximum value through the "bubble". Assuming uni- 
form conditions within the "bubble" for the magnetic field and 
the thermal electron density, we find the total amount of Faraday 
depth to be ~ 5 rad m"^, which agrees with the observed value. 
Finally a lower Faraday depth value for the nearby surround- 
ing "Fan" mediu m is wel l supported by previous observations 
( [Brouw & Spoelstra|1976| ). 

The chaotic and patchy appearance of observed features 
both in polarized intensity and polarization angle imply a large 
amount of depolarization and may be the result of the turbulent 
state of the probed magnetized medium. The power-law be- 
haviour of the power spectrum indicates a turbulent medium. In 
figure 14 we show power spectra of the bubble and background 
components; the spectral energy distributions are close to, but 
slightly flatter than, a Kolmogorov spectrum. Furthermore, the 
presence of an inertial range is apparent, which extends for 
more than a decad^^and flattens at modes k > 10^ . Such an 
increase in power on small scales agrees with the expected large 
amount of depolarization. 



6.5. The model in physical space 

These considerations lead us to a qualitative model in physical 
space as depicted in figure [15] The region nearest to the observer 



^ Due to the finite size of the field of ~ 12°, the smallest mode we 
can consider is kmin ~ 10^, while the largest mode kmax ~ 2.6 x 10^ is 



defined by the angular resolution. 
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Fig. 13. Grey scale images of polarization angle at selected Faraday depths at (p = [-7, -5, -3, -1] rad m ^ from top left to bottom 
right. The gray scale runs from +90° (black) to -90° (white). 



is bordered by a synchrotron emission layer that corresponds to 
the Local Bubble wall, about 70 pc away (the "curtain"). This 
region has a small Faraday depth and is therefore (partially) ion- 
ized and has a parallel magnetic field component towards the 
observer. The border region has enhanced synchrotron emission, 
which points to an increase in magnetic field perpendicular to 
the observer. 

Behind this region is an extended, magnetized region corre- 
sponding to the "ring", a volume of warm partially ionized and 
weakly emitting medium which contains the "bubble" compo- 
nent. The parallel component of the large-scale magnetic field in 
this direction is expected to point slightly away from us ( |Sun et 



[aL][2008 ), in agreement with these two components at negative 
Faraday depths. Since most of the diff'use synchrotron emission 
is thought to be emitted in the spiral arms and the spatial dis- 
tances we infer are within a few hundreds of parsec, we conclude 
that the overall detected PI emission originates primarily in the 
Local arm. 



6.6. Possible associations 



Circularly shaped, polarized objects without counterparts in 



other tracers have been observed before (e.g. Gray et al 



19981 



Uyaniker & Landecker 2002) . These observations differ in two 
ways from the previous ones. 



Firstly, we observe much more small-scale fluctuations in 
the polarization angle patterns. However, this can be caused by 
the difl'erent observing frequencies: Faraday efl'ects are much 
stronger at low frequencies, so that angle variations are expected 
t o be more eviden t. Secondly, [Gray et aL| ( |1998| ) and |Uyaniker &| 
Landecker ( 2002| ) explain their polarization structur es in terms 
of enhancement in electron density, while we (and Haverkom 
et al.|2003] ) argue that some of the structure must be due to en- 
hanced magnetic field. The importance of the magnetic eff'ects 
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Table 4. Main properties of selected stars within the "bubble". Data are taken from the SIMBAD database. 



Star ID 


Spectral type 


a 


6 


Distance 


Prop, motions 


Radial vel. 







(J2000) 


(J2000) 


[pc] 


[mas/yr] 


[km/s] 


WD 0314+64 


DA 


03 18 35.07 


+65 00 01.2 




110.0-169.0 




HIP 15520 


B2IV:e 


03 19 59.273 


+65 39 08.250 


234+24 


11.25-14.56 


-3.4 





Fig. 14. Power spectral density of polarized intensity of the 
''ring'' (stars) and ''bubble'' (circles) components as derived 
from a 6 degree inner square of the PI images integrated over 
the relevant Faraday depths. Also statistical errors and the 
Kolmogorov spectral slope are shown by \cr error bars and a 
black solid line, respectively. 



and their coupling with the particle density is supported by the 
ratio of thermal to magnetic pressure which is given by 



Fig. 15. Cartoon illustrating the proposed model for the ob- 
served distribution of Faraday depths. Grey areas indicate only 
emitting regions, and white blocks indicate regions both emitting 
and Faraday rotating. Intrinsic Faraday depth values are shown 
within each structure in the top panel, and observed Faraday 
depth profiles are built up in the lower panel for two representa- 
tive lines of sight: the dotted line crosses all the detected struc- 
tures, the dashed line misses the "bubble". Arrows in the top 
panel point out the magnitude and sign ofneB\\ 



ISmekTg 



(10) 



where B > B\\ is the total magnetic field strength, k the 
Boltzmann constant, and the electron temperature in K. 
For the background "ring" component, with Te = 8000 K, 
B\\,ring ~ 1-2 yuG, and ~ 3 x 10"^ cm"^, we find j3pr < 1.15, 
which indicates a non-negligible role played by the magnetic 
field with respect the thermal motion. For the bubble (with 
rie < ne,Bubbie and B\\ < B\iBubbie) we find a lower value Ppr < 6.6.2. H ii region 



sure could be explained by a very old remnant. However, no 
traces of shocks are found at other wavelengths, and the radial 
Faraday depth gradient within the bubble does not agree with 
the expected 0-profile for supernova remnants ( [Uyaniker et al. 
2002[ [Harvey-Smith et al.||2010| . Finally, a supernova remnant 
is a shell, while the bubble is thought to be a filled spherical 
structure. 



3.62 X 10-^7^^8000 K), indicating that magnetic eff'ects dom- 
inate thermal eff'ects over a wide range of temperatures (i.e. 
Te < 2.2 X 10"^ K); therefore the observed morphology can 
be shaped by the magnetic pressure of the surrounding ISM. 
Because of the low thermal electron density, the estimated Ppr 
value is far outside the conmionly observed range (Ppr > 1) for 
an Hii region (Harvey-S mith et al.|2 011) unless a strong gradi- 
ent of temperature with respect to the surrounding is considered. 
To address the origin of the bubble we now examine several pos- 
sible origins of this component: an old supernova remnant, an 
old planetary nebula, and a photo-ionized region. 

6.6.1 . Old supernova remnant 

The presence of a Faraday depth gradient across the bubble, the 
abundance of small scales and an evident polarized intensity pat- 
tern could indicate a supernova remnant. The absence of Stokes I 
emission, as well as the low ratio of thermal to magnetic pres- 



We consider the case of an Hii region in an unusually low- 
density medium and a possible association with the nearby B- 
type star HIP 15520 dws ~ (234 ± 24) pc ( [Van Leeuwen|2007| ). 
The distance, position and proper motion of this star (see Table 4 
and figure [16) lie within the constraints for the bubble. A link 
between this star and the Fan region was first suggested by 
^rschuur ( [1969 ), who pointed out the neutral gas depletion 
around the position of this star. He explained the disturbance in 
the H I by the motion of this star through the interstellar medium. 

However, an active H ii region is a strong source of Ha. We 
calculated the estimates Ha flux from the bubble if it is an H ii re- 
gion. The ionization flux (2o) for a a B2 type is Qq < 10"^^ ph s~^ 
( Ste rnberg et al.|2QQ3| , which can be used to constrain the den- 
sity in the Stromgren sphere as 



10-^ cm" 



1.2 pc 



?>li 



Go 



5 X 1049 photons s" 



1/2 



(11) 
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3'' 20"' m- 

RA (2000.0) 



Fig. 16. Relative position of selected stars (HIP 15520 is given 
by the cross, WD 0314+64 by the dot) with respect to the po- 
larized bubble here indicated by the circle. The offsets from the 
bubble centre (indicated by a plus) are ~ 0.23° and ~ 0.83°, 
respectively, while the bubble radius is ^ 1.15°. 



for it to be able to explain the bubble, it must have been much 
hotter and ionizing its surroundings when it was located near the 
centre of the current bubble. 

From the above reasoning, it is possible that the bubble is 
a fossil Stromgren sphere created out of an old planetary neb- 
ula and associated with a white dwarf star. Owing to the old age 
of the relic, interaction with the surrounding interstellar medium 
is expected, and the high symmetry of the observed morphol- 
ogy suggests that the surrounding medium is rather uniform. 
For a white dwarf at a distance of 200 pc, with a density of 
the bubble of < 0.18 cm"^, the required ionization flux Qq 
is 2o ^ 10"^"^ ph s"^ This corresponds to an emission measure 
EM < 0.34 pc cm"^, consistent with observations. 

From the Lyman continuum fluxes of| Panagia| (|1973[) and the 
white dwarf temperature and radius, we estimate the minimum 
temperature needed for sustaining the bubble: > 26000 K. 
The time needed for the white dwarf to cool from to the 
current value is ~ 10^ yr (roughly estimated from the cooling 
time in |Mestel| ( |1952| )), which does not agree with the age es- 
timated from its proper motion. Therefore, this particular white 
dwarf is probably not the cause of the bubble. However, other 
white dwarfs, at even greater angular distances from the centre 
of the bubble, may have created the bubble. 



At a distance of 234 pc and for an angular size of ~ 3°, we find 
a corresponding radius of 3.9 pc, a density of ~ 0.4 cm"^ and an 
emission measure of EM hip 15520 ~ 5.25 pc cm"^. This expected 
emission measure value is well above the WHAM sensitivity, 
thus it appears that a steady H 11 region is not likely to be respon- 
sible for the observed bubble. 



6.6.3. Old planetary nebula 

A more likely scenario would be provided by the final stage of 
the life of a star, forming a compact object. For a white dwarf 
star, some outer shells of materials have to be expelled in the 
surrounding ISM before the collapse by a weak shock. After 
this ISM sweeping stage, the hot star can ionize the surround- 
ing medium quickly, producing a photo-ionized region. Indeed 
evolved planetary nebula nuclei can have Lyman continuu m pho- 
ton luminosities comparable to those of early B stars ( Reynolds 
1986), enabling rapid ionization of the surroundings. The ob- 
servational work of T weedy & Kwitter (1996) , focussed on old 
planetary nebulae, indicates there is a class of nebulae that ap- 
pear to be filled, in agreement with the bubble. 

Indeed, the white dwarf star WD 0314-h648 (see Table 4 and 



figure [16]) on the edge of the bubble, moving outwards, is an 
obvious first candidate. Its distance and radial motion are not 
known. However, the white dwarf has a high proper motion, 
which can be due to a small distance to the observer, or due to 
a kick acceleration at its formation. If we assume that this white 
dwarf was born at the approximate centre of the bubble, we can 
obtain an estimate of its age as the time needed for the star travel- 
ling from the bubble centre. We find an age of > 15000 yr for the 
star WD 0314-h648 close to the edge of the bubble (for a proper 
motion and an estimated radial distance from the bubble centre 
of /?m* ~ 0.2 arcsec yr'^ ( [Lepine & Shara|2005| ) and r* ~ 0.83°, 
respectively). This age is much shorter than the recombination 
time, thus supporting the inferred high ionization degree. 

The eff'ective temperature for white dwarf WD 0314-1-648 is 
17500 K (Sion et al 1988). This temperature is so low that the 
white dwarf will not ionize its surroundings, in agreement with 
the lack of any bubble features around WD 0314-h648. However, 



6. 7. Comparison of the model at other wavelengths 

We compare our model to two radio observations of this field at 
higher frequencies, the 85 cm WSRT observatio n by|Haverkom 
|et al.| ( |2()(J3| and the PI survey at 21 cm by [WoUeben et al. 
( 2006 ). In our framework the bulk of the emission must build up 
within and beyond the foreground "ring" component, at a few 
hundred pars ec from the o bserver, in agreement with the con- 
clusion of WoUeben et al. ( 2006 ). Also, in agreement with their 
claim of multiple emitting components over a wide range of spa- 
tial scales, we suggest there are at least two emitting and Faraday 
thin components: the very nearby "curtain" component and the 
"ring" component; however, the "bubble" component is not seen 
in their polarized intensity map. We explain this lack of detection 
as due to the low synchrotron emissivity of the "bubble" feature 
at this frequency (the "bubble" does not manifest in emission) 
and to the low depolarizing impact it has on the background 
emission. For a maximum (ptubbie ~ 8 rad m"^ the expected mod- 
ulation of the polarization angle at 21 cm is ~ 21 degrees, but 
over the beam of the map the mean (pbubbie is < 4 rad m"^ provid- 
ing a variation in the polarization angle < 10 degre es. Interesting 
is also the comparison with the 85 cm data of Haverkorn et al. 
(|2003|). Inde ed, the standard RM values estimated by Haver korn 
et al.| ( [2003 ) are in good agreement with the Faraday depths of 
the "bubble" and "ring" components we present in this paper. 
Moreover, the spatial distribution of RM values matches the ob- 
served Faraday depths observed at 2 m. However, the PI map 
at 85 cm displays a ring-like feature instead of a filled one. This 
diff'erence is not likely to be explained by beam or depth depolar- 
ization (the resolution of PI maps is similar at 1 m and 2 m) but 
may be the consequence of the missing short spacing on Stokes 
Q and U maps at 85 cm, as also noted by Bemardi et al. ( 2009 ). 
Finally, we note the similarity of our model with the one pro- 
vided by [Bemardi et al. ( 2009| ), who consider the presence of 
three separated components: a synchrotron-emitting plus a sepa- 
rated Faraday -rotating components beyond the "bubble", which 
correspond to our spatially extended, synchrotron-emitting and 
Faraday -rotating "ring" component, which surrounds the "bub- 
ble". The main diff'erences from their model are the absence of 
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the very nearby "curtain" component and the position of the 
"bubble" with respect to the "ring" component. 

7. Summary and conclusion 

We used the RM- synthesis technique on WSRT radio polariza- 
tion data at ~ 2 m to reveal a complex Galactic synchrotron- 
emission foreground in a field in the Fan Region centred at 
(l,b) = (137°, 7°). We detected polarized signals with Faraday 
depth between -13 and -h5 rad m"^, which can be separated into 
three distinct components around -5, -2, and -h2 rad m"^, based 
on morphological consistency and coherence in Faraday depth 
space. For the first time, cross-correlation was applied to iden- 
tify and characterize these polarized structures in Faraday depth 
space. The low Faraday depth values of all components suggest 
nearby locations of the emission. 

The structure at [0,-h5] rad m"^ (called the "curtain") very 
likely corresponds to enhanced synchrotron emission due to 
compressed magnetic fields in the Local Bubble wall at around 
70 pc distance, in agreement with the location of the wall from 
optical polarized starlight data. The low but positive Faraday 
depth of this component suggests a parallel magnetic field com- 
ponent directed towards the observer. 

The second component centred at -1 rad m"^ (called the 
"ring" component) is located just behind the Local Bubble. It is 
expected to be spatially extended and Faraday-thick, and there- 
fore only partially detected in these data. 

The third component (the "bubble") is located within the ex- 
tended "ring". Its regular circular shape suggests relatively uni- 
form electron density and magnetic field structure, as is expected 
within the Fan region. The Faraday depth of the "bubble" in 
combination of its non-detection in Ha indicates that it cannot 
be created by density enhancement alone, but needs to have en- 
hanced magnetic field strength as well, contrar y to earli er de- 
tected regular polarization structures ( G ray et al.|1 998 ; Uy aniker| 
& Landecker||2Q02t |Uyaniker"et al., ,20031 ). Several possible as- 
sociations of the bubble with ISM objects we discussed, but the 
most likely explanation is the presence of a nearby (~ 200 pc) 
relic Stromgren sphere, associated with old unidentified white 
dwarf star and expanding in a low-density environment. 

Polarized intensity from both the "ring" and the "bubble" 
components show a power spectrum that has a power law, which 
is indicative of turbulence. Also, the discrete Faraday depth com- 
ponents suggest discrete, small-scale synchrotron emitting struc- 
tures in the ISM, as noted earlier by |Brentjens| ( |2011 ). 

This analysis shows that radio-polarimetry can detect mag- 
netized objects in the ISM that is not detectable in any other 
means. If fainter and Faraday thin structures lie along the line 
of sight toward the Fan region, this will be seen by deeper low- 
frequency observations. In the near future, the higher sensitivity 
and resolution of modern lo w-frequency arrays , such as the LOw 
Frequency ARray (LOFAR, |Heald et al.|2011|) may reveal these 
new components. The study of these spatially very extended 
and polarized foreground components will benefit from accurate 
Faraday rotation measurements (e.g. to separate the low cp "cur- 
tain" and "ring" components), as well as from the wide-field and 
accurate polarimetric imaging capabilities. A LOFAR observa- 
tion over its full bandwidth (about 200 MHz from LBA to HBA 
high) will improve the RM-synthesis maximum scale sensitivity 
in Faraday depth space, allowing to be detected also the eventu- 
ally Faraday thick components along the line of sight. The low- 
frequency arrays of next generation will sensitively probe weak 
magnetic fields with low Faraday depths, which may reveal more 
("bubble"-like) magnetized objects. 
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